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ABSTRACT 

Faraday rotation of radio source polarization provides a measure of the integrated magnetic field 
along the observational lines of sight. We compare a new, large sample of Faraday rotation measures 
(RMs) of polarized extragalactic sources with galaxy counts in Hercules and Perseus-Pisces, two 
nearby superclusters. We find that the average of RMs in these two supercluster areas are larger than 
in control areas in the same galactic latitude range. This is the first RM detection of magnetic fields 
that pervade a supercluster volume, in which case the fields are at least partially coherent over several 
megaparsecs. Even the most conservative interpretation of our observations, according to which 
Milky Way RM variations mimic the background supercluster galaxy overdensities, puts constraints 
on the IGM magneto-ionic "strength" in these two superclusters. We obtain an approximate typical 
upper limit on the field strength of about ~ 0.3 fiG Z / (500 kpc), when we combine our RM data 
with fiducial estimates of electron density from the environments of giant radio galaxies, and of the 
warm-hot intergalactic medium (WHIM). 

Subject headings: intergalactic medium — magnetic fields — plasmas — Faraday rotation - radio 
continuum:galaxies — large scale structure 



1. INTRODUCTION 

The detection of a magnetized intergalactic medium 
(IGM) at levels of a few /j,G within clusters of galaxies 
(Kim et al. 1990, Kim et al. 1991, Feretti et al. 1995, 
Clarke et al. 2001) raises the question of what IGM 
magnetic fields might be detectable in the wider inter- 
galactic medium, well beyond galaxy clusters themselves, 
but within galaxy overdense "filaments" of large scale 
structure (LSS). The recent, much improved definition of 
filament-like ga laxy distributions in the Sloa n Digital Sky 
Survey rSDSS ^Abazaiian et al.l27)Ml27)0l 120051 ^ . in ad- 
dition to the more nearby Tw o Micron All Sky Survey 
f2MASS1 (|2MASS Teaml2003T) and the secon d Center for 
Astro physics galaxy redshift survey fCfA2 UCfA2 Team! 
1998) at z^0.17 makes it easier to identify galaxy over- 
dense filaments. In the local Universe, to which the two 
latter surveys apply, the superclusters' larger angular 
extent encompasses an optimally large number of radio 
sources with measured Faraday rotations. 

This paper describes an attempt, using the best RM 
data so far available, to detect an RM signal due to an as- 
sociated intergalactic magneto-ionic medium in the Her- 
cules, Virgo, and Perseus-Pisces superclusters of galax- 
ies. As discussed below, different theoretical arguments 
and simulations cover a wide range of predictions, thus 
strongly motivating the task of obtaining possible con- 
straints from direct observations. 

2. ARGUMENTS FOR AND OBSERVATIONS IN SUPPORT 
OF A WIDESPREAD INTERGALACTIC 
MAGNETO-PLASMA 

2.1. Proposed early magnetic field generation and 
amplification in the course of LSS evolution 

Magnetic field seeding in the primeval intergalactic 

medium has been proposed by Harrison (1970), in which 

of which we describe in this paper. 
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chaotic motions of charged particles can seed an ini- 
tial magnetic fie ld by a Biermann Battery-like process 
l)Biermannl 1 1 9501) . More recent calculations by Gnedin, 
Ferrara & Zweibel (2000) proposed to explain how the 
breakout of ionization fronts from protogalaxies and their 
propagation through hypothesized dense neutral fila- 
ments in the evolving web of large scale structure can 
also generate, and amplify an intergalactic magnetic field 
before z ~ 5. Essentially (d\B\/dt) in the magnetic in- 
duction equation is positive, so that weak pre-stellar in- 
tergalactic fields in effect convert kinetic into magnetic 
energy. An IGM field during and after formation of the 
first galaxies is a natural consequence of early galaxy 
star-driven outflow, as modeled by Kronberg, Lesch, & 
Hopp (1999), especially earlier than z ~ 7. Analogous 
models by Furlanetto & Loeb (2001) for black hole/jet 
driven outflow have been shown to have a similar IGM 
magnetic field seeding effect at earlier epochs. 

Space- filling fields at early epochs can be considered as 
seed fields to incorporate into computational models of 
large scale structure (LSS) evolution. The intergalactic 
fields could be amplified in the course of gravitational col- 
lapse of baryonic and dark matter into the "filaments" of 
large scale structure, essentially converting gravitational 
infall energy into magnetic energy. The combination of 
shearing and shocks have been calculated by some au- 
thors (e.g., Ryu et al. 1998) to generate magnetic fields 
up to ~ 0.1 /iG in these intergalactic filaments. Up to 
now, there has been no observational verification of the 
actual field strengths in large scales in filaments. Indeed, 
independent LSS evolution simulations by other groups 
(e.g., Dolag et al. 2005) arrive at very low LSS IGM 
magnetic field strengths, too weak to be detectable in our 
present Faraday rotation data. Remaining uncertainties 
in the detailed physics of these difficult simulations un- 
derline the importance of direct observational tests, one 



2 



Xu ct al. 



2.2. Black hole seeded intergalactic magnetic fields and 

cosmic rays 

A very different source of IGM magnetization origi- 
nates in the reservoir of gravitational infall energy in 
galactic central black holes (e.g., Kronberg et al. 2001). 
The average comoving density of galactic black holes 
(BH) in the mature universe is estimated at pbh ~ 
2 x 10 5 M /Mpc 3 and they are widespread in the popu- 
lation of galaxies (e.g., Gcbhardt et al. 2000). When we 
combine this fact with the radio source evidence that cen- 
tral BH - driven jets and radio lobes convert a substantial 
fraction of the BH formation energy into magnetic fields 
that are transferred to the IGM via radio lobes (e.g., 
Kronberg et al. 2001), the BH-generated magnetic field 
strengths in galaxy walls and filaments is estimated to 
be of order 0.1 pG or more (Colgate and Li 2004). 

Independent calculations of the BH photon energy re- 
lease during the quasar epoch near z = 2 have been made 
by Soltan (1982), Small & Blandford (1992), and Chokhsi 
& Turner (1992). The efficiency of conversion of the pre- 
sumed central BH infall energy to photon energy released 
is remarkably high, a few tens of percent. Analyses of 
giant radio galaxies (GRGs) show that the conversion ef- 
ficiency of this BH infall energy into magnetic energy is 
comparably high (Kronberg et al. 2001, 2004), i.e. the 
magnetic energy, cosmic ray energy and photon energy 
released into the IGM are comparable. 

2.3. Estimates of a detectable LSS RM signal based on 

analyses of Mpc scale BH-powered radio sources 

The Faraday rotation through an intergalactic re- 
gion permeated by a magnetic field and thermal (non- 
relativistic) gas is 

RM = 0.81 • • — | ■ — rad m" 2 . (1) 

cm - ' 3 ^G pc 

Over a dimension L = 1 Mpc, containing a 1 /iG 
field and having = 10~ 6 cm~ 3 , the estimated RM 
- 1 rad m- 2 . The RM would be reduced by 0(N 1 / 2 ) 
where N is the number of field reversals over a line 
of sight through the supercluster. The outer lobes of 
GRGs (;>1 Mpc) contain at least partially aligned fields 
over at least a fraction of a megaparsec (e.g., Kronberg, 
Wielebinski & Graham 1986, Mack et al. 1997). Simi- 
lar evidence is provided by the peripheral regions of the 
synchrotron-emitting cosmic ray halos of some galaxy 
clusters (Clarke & EnBlin 2005). These facts suggest 
that some early component of large scale ordering ex- 
isted within typical supercluster volumes. 

Information on the plasma densities from direct ob- 
servation can be gleaned from detailed, multifrcquency 
observations of a few GRGs having dimensions up to ~ 
4 Mpc. Several such systems have established RM vari- 
ations that probe the "ambient" intergalactic medium 
within ~ 3 Mpc of the radio source and give estimates 
of n I t £\B IG \. Typically, n\% ~ lO^-lO" 5 cm" 3 for the 
ambient IGM on <~l-3 Mpc scales around the GRGs, and 
magnetic field strengths, based on the assumption of par- 
ticle/field energy equipartition are close to 1 /xG (e.g. 
Willis & Strom 1978, Strom & Willis 1980, Kronberg, 
Wielebinski and Graham 1986). Magnetic field strengths 
at these levels are independently consistent with values 
of <; 1 pG estimated by Kim et al. (1989) from very low 



levels of intergalactic synchrotron emission at 0.3 GHz 
that they detected in the IGM out to ~ 2 Mpc from the 
Coma cluster of galaxies. 

The above measurements lead us to estimate possi- 
ble RMi gm values from a supercluster environment that 
could be as high as ~20 rad m~ 2 in an optimistic 
scenario, for a typical supercluster pathlcngth scale of 
- 5 Mpc, n\% ~ 5 x 10~ 6 cm" 3 , and (\B IG \) ~ 1 /iG. 
It is more likely that RM contributions from LSS are 
typically less than ~ 20 rad m~ 2 , but some fraction of 
this RM could be detectable, depending on our ability to 
remove Galactic foreground RM, and this is what moti- 
vates us to attempt it. 

3. DATA AND ANALYSIS TOOLS 
3.1. Faraday Rotation Measure Data 

We have used an augmented version of the extra- 
galactic source Faraday rotation measures published by 
Simard-Normandin, Kronberg and Button (1981). The 
number of RMs that are suitable for our LSS filament 
probes is roughly double that in the Simard-Normandin, 
Kronberg & Button catalog of 1981. The accuracy of 
some of the RMs has also been improved in the present 
analysis by the addition of more radio polarization mea- 
surements. We omitted the very few RM sources that are 
in front of the local supcrclusters so as to have only back- 
ground RM "illuminators" of the supercluster medium. 

We define the supercluster boundaries using optical 
galaxy counts as explained further below. Because Fara- 
day rotations of radio galaxies and quasars at z ^ 0.5 are 
commonly affected by optically unseen intervening sys- 
tems and source-intrinsic evolutionary effects (Kronberg 
& Perry 1982, Welter, Perry & Kronberg 1984, Athreya 
et al. 1998) we have omitted the RMs of most sources 
that are at large redshift since they, as expected, show 
RM anomalies that are unrelated to the RM-weak, lo- 
cal universe intervenors that that we wish to probe. We 
also removed RMs of all sources lying within 7° of the 
Galactic plane because of the strong RM perturbations 
due to the Milky Way. Finally, since galaxy clusters are 
known to perturb the Faraday rotations of radio sources 
within and behind them (Kim ct al. 1990, 1991, Clarke 
et al. 2001), we have removed RM sources within 0.5° of 
known galaxy clusters. 

The increased numbers of Faraday RMs now available 
to us are sufficiently large that we have for the first time, 
relatively good sky coverage behind and beside the lo- 
cal superclusters. We used two RM subsets, (i) sources 
with known z < 0.5, and (ii) a comparable number of 
sources with unknown z, but from which we have re- 
jected a small number of anomalous values which could 
be due to high z effects or intermediate redshift clus- 
ters. The second sample is comparable in size to the first, 
and gives us substantially more sources at the expense 
of some possible high-z effect "pollution" . Consequently 
we combined the two subsets in all subsequent analysis, 
and have "cleaned" out statistically deviant high RMs. 

We distinguish between two kinds of RMs: the in- 
dividual measured source RMs and the smoothed RM 
(SRM). The latter are a combination of galactic fea- 
tures and, possibly, given the "expectation" arguments 
outlined above, degree-scale contributions from local- 
universe large scale structure. SRMs are obtained by 
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averaging individual RMs within a circle of a radius of 
15°, centered on the sources. The SRM smooths out 
small-scale RM fluctuations due to the radio sources, 
and should mainly contain information on magnetic fields 
from large angular scale structures. In order to dis- 
play the SRM variations more clearly and to compare 
with the superclusters visually, we Gaussian-smooth the 
SRMs to a 7-degree scale, to match the smoothing scale 
for the 2MASS data; A 7-degree smoothing is sufficient 
to clearly display the supercluster profiles. The upper 
panel of Fig. 1 shows a plot of the Gaussian-smoothed 
distribution of the radio source SRMs for the galactic 
sky that are relevant to the present analysis. The black 
points show the positions of the individual sources used 
in the SRM calculation. The colorbar gives the scale for 
our best estimate of the SRM over the sky. In order to 
more clearly illustrate the contrast of the lower SRMs, 
we saturated the SRM color scale at —65 rad m~ 2 and 
+65 rad m -2 in blue or red respectively. The larger 
SRMs extended to ~ 1 100 1 rad m~ 2 . The southern sky 
below 5 =-25° is excluded because of the lower angular 
density of RMs, hence less well-determined SRMs. We 
also exclude the area near the galactic equator, i. e., 
| 6| < 7°, due to the strong galactic con tribution to SRM 
( Sim ard-Normandin fc Kronberdll980|) . 

3.2. The CfA2 and 2MASS Galaxy Surveys 

We have utilized two sets of galaxy survey data for our 
analysis, the second Center for Astrophysics survey and 
the Two Micron All Sky Survey. 

The second Center for Astrophysics (CfA2) Rcdshift 
Survey covers the whole northern celestial hemisphere, 
and has measured redshifts of about 18,000 bright galax- 
ies with apparent blue magnitude mf, less than 15.5. This 
survey captures the three-dimensional properties of the 
large-scale structures in the nearby universe (z < 0.05), 
whi ch is very useful fo r RM-LSS cross-correlation analy- 
sis (|CfA2 Teamlll998j) . 

The Two Micron All Sky Survey (2MASS) has cov- 
ered 99.998% of the sky observed from both northern 
and southern 2MASS observation sites l)2MASS Teaml 
2003). The 2MASS has three infrared cameras in the J 
(1.25 £tm), H (1.65 /an), and Ks (2.17 /zm) bands. Near- 
infrared observations are useful for lines of sight near the 
Galactic Plane since the extinction due to the interstellar 
medium at near-infrared wavelengths is only about one- 
tenth of that at optical wavelengths. The all-sky data 
release from 2MASS includes 470,992,970 point sources 
and 1,647,599 extended sources or galaxies in the nearby 
universe. In order to keep the data uniform, we only 
use those extended sources with apparent Ks-band mag- 
nitude brighter than 14.5 (i.e., m,Ks < 14.5). This still 
leaves 1,312,603 extended sources for our analysis. 

The 2MASS contains only two-dimensional informa- 
tion of LSS, while the CfA2 contains three-dimensional 
information. Therefore, as described in Subsections 14.11 
and 14.21 in detail, we use different tools to analyze the 
data from the two surveys. For 2MASS, we compare the 
column density with RMs, calculated using the HEALPix 
algorithm (section 4.1), whereas for CfA2, we use the 
Voronoi tessellation method, described next. 

3.3. 3D Voronoi Tessellation 
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Voronoi Tessellation (VT) is a robust tool for cal- 
culating number densities in two-, three- or higher- 
dimensional spaces. Here, we use a very efficient VT al- 
gorithm, QHULL (Barber, Dobkin & Huhdanpaa 1996), 
to calculate the galaxy number density in the CfA2 cat- 
alog. The procedure is summarized below. 

First, using QHULL, we construct the entire volume 
of the CfA2 catalog with a number of Voronoi Diagrams 
(VDs). Each VD corresponds to a galaxy in this cata- 
log. The VD of a given galaxy is a convex hull bounded 
by a series of planes that bisect the segment between 
this galaxy and its neighbors. In consequence, any query 
point inside this VD is closer to this galaxy than any 
other galaxy. The "proprietary volume" of a galaxy is 
that of its VD (|Icke & van de Wevgaertlll987j) . The lo- 
cal galaxy number density is the inverse of its volume, 
i.e., pi = 1/Voh, where pi and VoL are the number den- 
sity and volume for the i t/l -galaxy. The overdensity of a 
galaxy is defined as 

h = (pi ~ Pi)/pi, ( 2 ) 

where pi is splined from p(z), which is the mean number 
density at a given redshift bin at z. The rotation mea- 
sure should relate to both the galaxy number density 
and path length of a radio ray passing through a por- 
tion of large-scale structure volume. Due to the survey 
selection function, the mean number density p depends 
strongly on galaxy redshift. In order to eliminate the 
effect of this selection function, we employ an estimator 
of the galaxy overdensity by effectively normalizing out 
the survey selection function. The overdensity describes 
the relative fluctuations of galaxy density, therefore in 
principle it should not depend on redshift. At this point 
we are in a position to use the overdensity of galaxies for 
our correlation analysis. 

3.4. Pathlength and weighted pathlength 

Aside from the overdensity, another important param- 
eter is the pathlength of radio sources through the super- 
clusters. First, we need to ascertain whether the line of 
sight to a radio source passes through superclusters. Sec- 
ond, we need to calculate the effective pathlength. Since 
we have tessellated the superclusters into a large num- 
ber of VDs, the pathlength to a radio source is the sum 
of lengths inside the VDs that radiation from the source 
passes through, i.e., 

plj = ^2hj, (3) 

i 

where plj is the pathlength, and Uj is the length inside 
the i th VD which radiation from radio source j passes 
through. We also introduce a related parameter, the 
weighted pathlength, wplj, which is the sum of the prod- 
uct of overdensity and length inside VDs: 

i 

where 5ij is the overdensity of the i th VD which the radio 
source j passes through. As defined in Eqn. [21 overden- 
sity is a linear function of the real galaxy number den- 
sity, Hence, the product of overdensity and pathlength is 
a linear function of galaxy column density at the specific 
VD. 
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Fig. 1. — (a) All-sky plot of the extragalactic source SRM distribution. The square black points show the positions of individual SRM 
sources. The color map shows SRM variations Gaussian-smoothed to 7 degrees in an Aitoff projection of galactic coordinates. The color 
range from blue to red signifies a range of SRM values from —65 rad m~~ 2 to 65 rad m -2 . The excluded areas (see text) are for |6| <7° 
and 8 <-25° (b) The 2MASS galaxy column density Gaussian-smoothed to 7 degrees in the same coordinates as in (a). The Hercules and 
Perseus-Pisces superclusters lie in the high density parts labeled I and IV, respectively. Labels II, III and V indicate the RM comparison 
zones (see subscction lij.il and table HI . 



4. RESULTS AND DISCUSSION 

In order to correlate the RMs of radio sources with 
galaxy overdensity in nearby superclusters, it is crucial 
to find a way to identify those superclusters or local fila- 
ments of large-scale structure. We have developed a fast 
algorithm to identify superclusters after tessellating all 
galaxies in the Cf A2 catalog . The algorithm will be de- 
scribed elsewhere ijXul 12005). The basic idea is to find 
structures (or continuous Voronoi Diagrams) with over- 
densities in a certain range. For CfA2, we found that an 
overdensity larger than 2 sets a very reasonable thresh- 



old to identify superclusters. Using this algorithm, we 
successfully identified and defined 9 nearby superclustcr 
candidates. 

Out of 9 supercluster candidates that we initially iden- 
tified in the CfA2 catalog, only the Perseus-Pisces, Her- 
cules, and Virgo superclusters give a sufficiently large 
angular coverage in a favorable part of the galactic sky, 
with a sufficiently good match of galaxy column density 
and angular RM coverage on the sky. This is needed to 
ensure that we have a sufficient number of radio sources 
in those sky areas to statistically analyze the data. The 
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Virgo supercluster is the least favorable of these three, 
being inconveniently close to us, hence distributed over 
a very large area of sky, so that galactic foreground RM 
effects are more difficult to remove. Further, the low 
angular number density of Virgo supercluster galaxies 
is relatively more contaminated by background galaxies. 
Therefore, we compare the RMs with 2MASS data only 
for the Perseus-Pisces and Hercules superclusters, al- 
though we did attempt the VD method for Virgo (Fig. 2). 

In this analysis it is important to recall that "fore- 
ground" SRM features due to the Milky Way probably 
exist at some level at the higher galactic latitudes where 
our superclusters lie. Our search for an RM signal from 
local supercluster fields can therefore only rely on de- 
tailed correlations between zones of galaxy overdensity 
and of the SRM. The (reasonable) premise here is that 
Milky Way-induced SRMs are not correlated with extra- 
galactic magnetized regions. Even if galactic + inter- 
galactic superpositions exist at weak RM levels we can, 
at minimum, set the first useful upper limits to SRM 
contributions from supercluster zones. Such upper lim- 
its are interestingly close to the predictions of theoretical 
calculations, as discussed in Section 3. 

There is an intriguing juxtaposition of a Milky Way 
and extragalactic RM contribution in Region A, a well- 
known anomaly in the RM sky (Simard-Normandin & 
Kronberg 1980). Region A, the huge negative RM zone 
located in 50° < I sC 150° and b ^ -40° [see the large 
deep blue area in Fig. 1(a)], overlaps in part with the 
Perseus-Pisces supercluster, but has a greater overall 
angular extent. As we discuss below, there is a weak, 
but positive correlation (Table 1) between the SRM and 
galaxy overdensity in this zone. If this correlation is not 
being (perversely) masked by galactic effects, it would in- 
dicate that ordered magnetic fields in the Perseus-Pisces 
supercluster maintain the same prevailing sense on mega- 
parsec scales. Unfortunately firm conclusions on this 
question cannot yet be drawn because of possible "in- 
terference" by the Milky Way foreground RM. Alterna- 
tively stated, the P-P supercluster filament "rides" on 
the large negative RM zone of region A. One possible 
interpretation is that the P-P supercluster contribution 
to the SRM (Fig. 1) contains a positive RM component 
that is not strong enough to overwhelm the galactic RM 
contribution in this same area. 

The RM region A does not appear to correlate with 
other known large scale features of the galactic sky, 
such as the HI column density, diffuse galactic syn- 
chrotron emission, gamma rays, etc. The galactic 
feature that does generally coincide with region A is 
the diffus e 45/zm radiation, due presuma bly to galactic 
dust l|Schlegel. Finkbeiner k Davis! H998). It is further 
salutary that the RMs in Region A are predominantly 
negative, which conforms to the prevailing direction of 
the interstellar field in the Sun's vicinity. However this 
fact by itself need not necessarily refute the existence of 
a mixed galactic/extragalactic origin for the anomalous 
Region A. 

4.1. Correlation between RMs and 2MASS Data 

2MASS is an all-sky photometric survey with three in- 
frared cameras. An easy way to calculate the 2MASS 
galaxy number density is to pixelize the all-sky data 
appropriately. Pixels corresponding to the location of 



large-scale structures should have a higher number den- 
sity of galaxies. The public software package Hierar- 
chical Equal Area isoLatitude Pixelization ( HEALPix) 
of the sphere l|G6rski. Hivon. &: Wandeltlll999j) is well- 
suited to this purpose. We pixelize the 2MASS data into 
12 x 128 2 cells. Since each pixel has the same area, we 
need concern ourselves only with the number of galaxies 
in each pixel. In order to obtain a continuous density 
field, we smooth the galaxy column density {N g ) onto a 
1° scale to eliminate discontinuities on very small scales. 
After smoothing, the mean galaxy number in each pixel is 
around 7. Therefore, we classify the number density into 
three ranges: an underdense range (N g < 7), one sigma 
higher than the mean density range (7 ^ N g ^ 10), and 
an overdensity range (N g > 10). 

To illustrate the nearby superclusters more clearly, we 
Gaussian-smooth N g on a scale of 7 degrees as shown 
in the lower panel of Fig. 1. The Hercules and Perseus- 
Pisces superclusters constitute the high density parts la- 
beled I and IV, respectively. Labels II and III indicate 
the (low galaxy density) RM comparison zones for region 
I in the same Galactic latitude ranges (30° ^ b ^ 70°); 
Label V indicates the comparison zone of region IV 
(-40° ^ b < -15°). We do not use region VI since 
the RM data there are less complete. 

As estimated in Section 12.31 RM values from a super- 
cluster environment might be as high as ^20 rad m~ 2 . 
If we assume that N g is proportional to the pathlength 
of radio sources through large-scale structures, and that 
is constant, we would expect positive correlations 
between RMs and N g . Due to the unknown complexity 
of alignment and reversal scale of the IGM magnetic field, 
it is hard to find a direct correlation between the abso- 
lute RMs and galaxy counts in the corresponding pixels. 
However, we may expect that the mean values of absolute 
RMs in a supercluster area should be larger than in other 
areas without superclusters. Due to the magnetic field 
of the Galaxy, however, observed SRMs strongly depend 
on th e galactic latitude ijSimard-Normandin &: Kronberd 
1980). We therefore analyze the RM data in two corre- 
sponding spherical rings (see Tableland Fig. 1). Tabled 
lists mean values of absolute SRM of the various areas 
classified above. We have filtered out zones of anoma- 
lously high local galaxy density, so as to avoid contami- 
nation by e.g. galaxy clusters. 

For the Hercules supercluster zone we find that over 
its entire area the mean SRM is consistently larger than 
in the two comparison zones. It is close to, but does 
not coincide with the North Galactic Spur (Haslam et 
al. 1982). This systematic difference also holds for each 
subrange of galaxy density (Table 1). However, whereas 
there is a global correlation of enhanced SRM with the 
overall boundary of the Hercules supercluster, the de- 
tailed locations at 7° resolution of maxima in smoothed 
galaxy column density with those of the SRM [compar- 
ing Figs. 1(a) and (b)] do not correspond exactly. This 
is also reflected in the plot of weighted pathlengths (Sec- 
tion ^. 4|) against RM in Fig. 2, in which we used the VD 
method and the CfA2 data fSection 14. 2|) . We can in- 
terpret the combination of these two results to indicate 
that an RM component is generated in the Hercules su- 
percluster, but if so, most ordered magnetic field regions 
avoid the zones of highest galaxy density. This conclu- 
sion must remain tentative until we better understand 
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TABLE 1 

Comparison of the mean values of absolute SRM in areas I, II, and I+II+III, and in areas IV and V. 



Sky Area 


Galaxy Density Range SRM mean 
[rad m -2 ] 




all 


14.3 ± 2.0 


rlGrculcs J\yg3i 1, oU sources 


Ng<7. 


14.8 


rn° < / < 7fi° sn° < k < 70°) 


7. < Ng < 10. 


15.9 




Ng > 10. 


12.9 




all 


10.0 ± 1.8 


Hercules Comparison Area, II 


Ng < 7. 


8.7 


(70° ^ I ^ 120° or I ^ 330°. 30° ^ b ^ 70°), 30 sources 


N > 7 


10.9 




all 


10.2 + 0.8 


Hercules Comparison Zone, Area I+II+III, 154 sources 


Ng <7. 


11.4 


(0° < I < 360°, 30° O ^ 70°) 


7. sC Ng <c 10. 


8.7 




Ng > 10. 


10.5 




all 


56.8 + 8.6 


Perseus-Pisces Area IV, 45 sources 


Ng < 7. 


65.7 


(100° ^ I ^ 150°, -40° ^ b ^ -15°) 


7. < Ng <C 10. 


53.5 




Ng > 10. 


55.5 




all 


41.0 + 6.3 


Perseus-Pisces Comparison Area V, 42 sources 


Ng < 7. 


49.4 


(150° ^ I ^ 180° or I ^ 100°, -40° ^ b ^ -15°) 


Ng>7. 


32.6 



Note. — As illustrated in Fig. 1, areas I and IV are high galaxy density areas where the Hercules and Perseus-Pisces superclusters are 
located respectively, and areas II and V are low galaxy density areas adjacent to these two superclusters. 



the galactic foreground RM (which may also play some 
role, as suggested by the galactic pulsar RMs) . A future 
higher density of pulsar RMs is most crucially needed 
to clarify the decomposition of galactic and superclustcr 
RM contributions in this high latitude zone. 

In the case of the Perseus-Pisces supercluster, the 
SRMs within the overall supercluster boundary are sim- 
ilarly larger, again for all galaxy density regimes (Ta- 
ble 1). In contrast to the Hercules supercluster, there 
is also a suggestion of a positive correlation on smaller 
scales between SRM and galaxy density. This is seen 
both in Fig. 1, and in the plots (Fig. 2) of wpl vs. SRM. 

For the Virgo supercluster, we find no evidence for a 
correspondence on any smoothed scale (Figs. 1 and 2). 

4.2. Correlation between RMs and CfA2 data 

The CfA2 galaxy redshift survey covers the northern 
celestial sky. Using the Voronoi Tess e llation technique 
mentioned in Section 13.31 and in iXul l)2005[K we iden- 
tify galaxy superclusters in the CfA2 catalog. Then, 
from Eqns. 13141 we calculate the real pathlength pli and 
weighted pathlength wpli of the i-th radio source passing 
through the supercluster. Although wpli has the dimen- 
sion of length, it is a linear function of galaxy number col- 
umn density, which we assume is proportional to the col- 
umn density of thermal electrons in the IGM. The RM is 
proportional to the electron column density and the line 
of sight component of the magnetic field (Cf. Eqn. [TJ. 
Therefore, we expect a positive correlation between RMs 
and wph . Fig. 2 shows plots of weighted pathlength wpl 
vs. source RM and SRM respectively. The left three 
panels, from top to bottom, show the relation of wpl and 
source RM in Hercules, Perseus-Pisces, and Virgo super- 
clusters respectively, and the right three panels show the 
relation oiwpl and SRM in those three superclusters. We 
find that the wpl is positively correlated to the fluctua- 



tions of absolute SRM in the Perseus-Pisces supercluster, 
but that there are no such correlations in the other two 
superclusters (see above). 

4.3. RM data vs. RM expectations in the Hercules 
supercluster 

In Section l2~3l it was roughly estimated that the RM 
in a supercluster may reach ~ 20 rad m -2 , scaled to 
pi ~ 5 Mpc, n{% ~ 5 x 1CT 6 cm" 3 , and {\B IG \) - 1 fiG, 
by using the simple model described in Eqn.^and the re- 
lated text in Sect ion l2~3l RMs estimated from this model 
are designated as MRMs. Comparing these MRMs with 
the measured SRMs, we put constraints on the IGM mag- 
netic fields and electron density. The detailed procedure 
is described as follows. First, we use the simple model 
to estimate MRMs of all radio sources in the Hercules 
supercluster for different magnetic field reversal scales 
by varying n\® and \Bjg\. Second, we calculate the x 2 
of each given pair of and |-Bjg|, for a given i?-field 
reversal scale. Specifically, 

X 2 = Y^iMRM* - SRMi) 2 /a 2 SRMi , (5) 
i=i 

where N s is the total number of radio sources pass- 
ing through the Hercules supercluster. Third, as usual, 
the corresponding normalized likelihood is defined as 
C = cxp (-(x 2 - xLJ/ 2 ); where the Xmin is the min- 
imum x 2 °f a certain pair of and \Bjg\. Finally, 
we may derive the la, 2a, and 3a contours from likeli- 
hood values equal to 0.607, 0.135, and 0.011 respectively 
as shown in Fig. 3. From the calculated likelihood, we 
obtain constraints on the magnetic field and the ther- 
mal electron density. The four panels in Fig. 3 show 
the constraints on |Bjg| and for magnetic field re- 
versal scales of 200, 400, 600, and 800 kpc respectively. 
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Fig. 2. — Plots of weighted pathlength wpl vs. source RM and SRM respectively, from our VD-based analysis of the CfA2 catalog. 
The left three panels, from top to bottom, show the relation of wpl and source RM in Hercules, Perseus-Pisces, and Virgo superclusters 
respectively, and the right three panels show the relation of wpl and SRM in those three superclusters. 



The constraints have the expected dependence on the re- 
versal scales, the larger reversal scales yielding a lower- 
magnetic-field constraint with a given electron density. 
This follows straightforwardly from Eqn. ^ Fig. 3 il- 
lustrates how the B-field, reversal scale, and are 
strongly degenerate. Independent n\2 measurements are 
needed to help break this degeneracy. For instance, the 
most recent warm-hot IGM (WHIM) obs ervations imply 
a hydro gen density of ~ 5 x 10~ 6 cm -3 (|Nicastro et alJ 
2005©. If we combine this global estimate with our re- 
sults here, we obtain i?-field constraints of 0.4 ± 0.2 /j,G 
(2a) and 0.3±0.1 /iG (2a) with reversal scales of 400 kpc 
and 800 kpc respectively. Here, we use the SRMs as 
distinct from the RMs generated by the medium in the 
filaments. Since the SRMs may contain a Milky Way 
component, the absolute values of SRMs are statistically 
larger than the RMs generated in the filaments. In this 
sense, the constraints shown in Fig. 3 should be consid- 
ered as upper limits for the magnetic fields in filaments. 

5. SUMMARY AND CONCLUSION 

We have performed the first search for intergalactic 
magnetic fields in three galaxy-overdense zones of the 
nearby universe: the Virgo, Hercules, and the Perseus- 
Pisces superclusters of galaxies. We have compared our 
results with simplified model expectations of IGM mag- 
netic fields and intergalactic thermal gas, the latter being 
represented by the WHIM. Algorithms were developed to 
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Fig. 3. — Constraints on magnetic field and thermal electron den- 
sity in the Hercules supercluster for magnetic field reversal scales of 
200, 400, 600, and 800 kpc respectively. The contours in each plot 
show the la, the 2tr, and the 3c consistency contours respectively. 



calculate zones of galaxy overdensity in these three ar- 
eas using the 2MASS and CfA2 galaxy survey data, and 
combined with new Faraday rotation measures. 

Our results are summarized as follows: 

(1) The large angular extent of the Virgo supercluster 
combined with its relatively low galaxy surface density 
make any Faraday RM contribution difficult to distin- 
guish from small variations in the foreground RM of our 
Galaxy. We obtain a null result with the best available 
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data on RM and galaxy counts, which have inadequate 
signal-to-noise for the purpose. 

(2) For the Hercules supercluster, the Faraday RMs are 
enhanced overall within the supercluster boundary on the 
sky. However within the boundaries, the enhanced RMs 
do not generally occur precisely at zones of maximum 
supercluster pathlength. This may reflect the existence 
of large scale coherence of extragalactic magnetic fields 
within the Hercules Supercluster, but outside of the most 
galaxy-overdense zones within the supercluster volume. 

(3) The Perseus-Pisces supercluster, which coincides 
with a larger "region A" , a zone of anomalously large 
RMs off the Galactic plane (Simard-Normandin and Kro- 
nberg 1980), shows the clearest detailed correlation with 
galaxy column density and supercluster volume path- 
length for background Faraday- rotated radio sources. 
This correlation makes it likely, but not absolutely cer- 
tain that some of the enhanced RMs are actually gen- 
erated in the IGM of the supercluster over and above a 
Milky Way contribution. As with the Hercules super- 
cluster, future larger numbers of pulsar RMs are needed 
to disentangle the relative RM contributions from the 
Milky Way and the ICM of the Perseus-Pisces superclus- 
ter/filament zone. 

(4) Comparison with simple models of "expectation" 
IGM fields gives a result consistent with typical IGM 
fields in the supercluster IGM of order 0.3 L^ 00 ^ pc (J-G at 
a 2a significance level. 

Future advances in probing the Faraday Depths ijBurnl 
119661 iBrentiens fc deBruvnl 12005^ of supercluster envi- 



ronments will require (1) much denser sky coverage 
of discrete sou r ce Fa raday rotation measurements as 
iGaensler et alJ l)2005j) have recently achieved for the 
Magellanic clouds, and (2) a better specification of the 
galactic magneto-ionic environment at off-plane direc- 
tions using deeper surveys of pulsar RMs and dispersion 
measures, and multi-frequency all-sky radio polarization 
surveys. In addition, we require (3) fainter spectroscopic 
X-ray detections of intergalactic gas transitions to bet- 
ter specify the temperature and density of the WHIM. 
The first two will be easily accomplishable with the next 
radio telescope ge neration such as the prop osed "Square 
Kilometer Array" ((Carilli fc Rawlingsl2004D . At present, 
the SDSS and other low redshift optical galaxy surveys 
are the most complete portion of the entire suite of ob- 
servations required to better characterize the physics of 
the IGM. 
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